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We experimentally resolve several weakly coupled nuclear spins in diamond using a series of novelly 
designed dynamical decoupling controls. Some nuclear spin signals, hidden by decoherence under 
ordinary dynamical decoupling controls, are shifted forward in time domain to the coherence time 
range and thus rescued from the fate of being submerged by the noisy spin bath. In this way, more 
and remoter single nuclear spins are resolved. Additionally, the field of detection can be continuously 
tuned on sub-nanoscale. This method extends the capacity of nanoscale magnetometry and may be 
applicable in other systems for high-resolution noise spectroscopy. 


I. INTRODUCTION 

Detection of single nuclear spins is an outstanding is¬ 
sue in magnetic resonance spectroscopy and imaging. It 
would be beneficial to molecular structure analysis, and 
will have a far-reaching impact on chemistry, biology, and 
medicine [1, 2]. Besides, nuclear spins are valuable as 
quantum registers for their long coherence times, partic¬ 
ularly in relation to diamond defects [3-8] . However, it is 
challenging to detect single nuclear spins owing to their 
weak magnetic moments and the noisy environment of 
solids. Recently, single nuclear spin detection has been 
studied in various systems, including electrical transport 
measurements of ^®®Tb in a single-molecule magnet [9] 
and all-electrical detection of in silicon [10] . Addi¬ 
tionally, using an individual negatively charged nitrogen- 
vacancy (NV) center in diamond, numerous protons in 
organic samples [II, 12], four ^®Si in silica [13], and a 
single surface proton [14] have been detected. 

In this work, the single electron spin of NV center in 
diamond is used as a magnetic sensor [Fig. 1(a)] [15-17]. 
Strongly coupled nuclear spins (with coupling strength 
exceeding the electron spin dephasing rate ~ 1 /T ^) can 
be well resolved through the electron spin level split¬ 
ting [18-20]. To resolve weakly coupled nuclear spins 
(with coupling strength comparable or less than the elec¬ 
tron spin dephasing rate), the environmental noise has 
to be suppressed and the sensor sensitivity for nuclear 
spin signals has to be enhanced. Dynamical decoupling 
(DD) methods [21], such as the widely used Carr-Purcell- 
Meiboom-Gill (CPMG) sequences and XY sequences, are 
the most commonly used tools at present. By periodi¬ 
cally flipping central electron spins, weakly coupled sin¬ 
gle ^^G nuclear spins [22-24] and clusters [25, 26] around 
NV centers have been detected. In this work, we ap¬ 
ply novel DD controls [27] which have advantages in re¬ 
solving remote single nuclear spins. Applying the well- 
designed DD control sequences, some nuclear spin signals 
greatly diminished by decoherence are shifted forward in 


time domain to the coherence range, and thus rescued 
from the fate of being submerged by the noisy spin bath. 
With this benefit, two remoter nuclear spins, nearly silent 
under CPMG controls, pronounce their existence under 
appropriately designed DD controls in our experiment. 
Besides, the resolutions for characterizing some single nu¬ 
clear spins can be improved using this method. 


II. THEORY 


The dynamics of nuclear spins around an NV center re¬ 
sults in the modulation of the electron spin coherence and 
causes the decoherence effect [28, 29]. This phenomenon 
can be viewed in a semiclassical picture, where the co¬ 
herence L(t) of a two-level system in a time-dependent 
noise field is approximately expressed as [30] 


L{t) ~ exp 
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with t being the total evolution time, S{uj) the noise spec¬ 
trum [Fig. 1(d)], and F{ujt) the filter function [Figs. 2(b- 
d)] associated to the applied control pulse sequence. DD 
controls make the electron spin selectively sensitive to sig¬ 
nals at specific frequencies and suppress unwanted back¬ 
ground noise according to the features of F{ujt). Nu¬ 
clear spins around the NV center induce discrete peaks 
on S{oj) as Fig. 1(d) shows [31], and thus can be ex¬ 
perimentally identified via the characteristic dips on the 
electron spin coherence [22, 27]. Note that the semiclas¬ 
sical treatment here is qualitative and heuristic, and the 
quantitative and rigorous calculation should be based on 
quantum mechanics. 

The filter function F{u!t) encapsulates the influence 
of DD control on decoherence, so it relates to the posi¬ 
tions and magnitudes of nuclear spin signals in the time 
domain. The filter function F{ujt) of a typical DD con¬ 
trol has a series of peaks, which center at ujtj{2Tr) = 





2 


(a) 





(b) 


\b 

f t 

MW! D 


\m = 0> 


(c) 



a natural but rough way to alter F(ujt). In fact, we can 
design DD control sequences by changing pulse distribu¬ 
tions in a systematic way to carry out desirable tuning of 
F{u)t). By doing so, the signals of weakly coupled nuclear 
spins with different distances and orientations could be 
selectively optimized by such method. 

A CPMG-n sequence is n repetitions of a | - tt - | unit. 
In this paper, we expand the repetition unit to contain 
three tt pulses. The three tt pulses within a repetition 
unit is symmetric but non-uniformly distributed, i.e., the 
second tt pulse is located at the center of the unit, and the 
first and third tt pulses are separated from the second one 
by 3rr with the parameter r (0 < r < 0.5) characterizing 
the relative positions of the pulses. In short, the new 
sequence is n/3 repetitions of a - tt - 3rr - tt - 

3rr - TT - unit [Fig. 2(a)]. In this case, the filter 

function depends on both n and r, and is expressed as 
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FIG. 1. (color online). Sketch of the experimental system and 
methods, (a) The experimental system consists of an NV cen¬ 
ter and several weakly coupled nuclear spins, (b) Ground 
state energy levels of an negatively charged NV center in an 
external static magnetic held B, which is 27 G and parallel to 
the NV axis. The zero held splitting D is 2870 MHz. We en¬ 
code the quantum transition |ms = 0) ■f->' |ms = 1) as a probe 
qubit and manipulate it by resonant microwave pulses, (c) 
The pulse sequences contain two green laser pulses for initial¬ 
ization and readout of the electron spin state with microwave 
pulses in between to control the spin, (d) Noise spectrum of 
the nuclear spin bath. The hve red circles correspond 
to the hve single nuclear spins resolved in the experiments. 
The blue circles correspond to other bath spins, which are 
randomly generated by a simulating program considering the 
experimental conditions. 


where n is divisible by 6. The dominant peaks of this 
function center at cot/ (27r) = (2fc— l)n/6 (fc = 1, 2, 3, • • •) 
for r 7 ^ 1/3. In particular, the first one, which is espe¬ 
cially responsible for nuclear spin detection, moves for¬ 
ward and comes out at a;t/(27r) = n/6. The height of 
these peaks, corresponding to the amplification of the 
signals, depends on the value of r. The height of the 
fc-th dominant peak [Fig. 2(f)] is 

hn{r) = {1 - 2cos [{2k - l)7rr]}^ . (3) 

Such a degree of freedom provides advantages in nuclear 
spin detection. In the following we present the results for 
n = 30. 


{2k — l)n/2 {k = I, 2, 3, • • •) for an n-pulse CPMG con¬ 
trol (referred to as GPMG-n). As the total evolution 
time t increases, the peaks sweep over the noise spec¬ 
trum from high frequencies to low frequencies. When a 
peak of F{u}t) coincides with a discrete peak in the noise 
spectrum, a coherence dip appears as qualitatively indi¬ 
cated by Eq. (I). The depths of dips depend on the noise 
features of target nuclear spins and the peak heights of 
F{u}t). For a sequence with n flipping pulses, the peak 
height scales as v?. In the case of small n, the peak is 
so short that some dips will be quite shallow and can 
hardly be identified. In the opposite situation (large n), 
the peak is unnecessarily tall. The dips may be oscil¬ 
lating and malformed, and the signals from individual 
nuclear spins will be blurred. So the peak heights of the 
filter functions should be properly adjusted. However, 
conventionally used periodic DD controls have limited 
freedom in tuning F{ujt). Increasing the pulse number is 


III. EXPERIMENTS 

We measured the spin coherence of an NV sensor em¬ 
bedded in a bulk diamond with nitrogen impurity of low 
concentration (< 5 ppb) and ^^G isotope of natural abun¬ 
dance (1.1%). Several weakly coupled ^^G nuclear spins 
[Fig. 1(a)] impose a.c. magnetic noise on the NV sensor 
with characteristic frequencies that depend on the exter¬ 
nal static magnetic field and hyperfine couplings. Typ¬ 
ically, the a.c. noise from different nuclear spins induce 
different modulations on coherence curves of the NV sen¬ 
sor. The sensor coherence under GPMG-30 (i.e., r = 1/3) 
control is illustrated in Fig. 2(g). The dips in the zones 
labeled by I, II, and III originate from three distinct nu¬ 
clear spins, respectively. We notice that, with GPMG-30 
(i.e., r = 1/3) control, the oscillating patterns in / and II 
blur the signals from individual nuclear spins. By apply¬ 
ing the designed DD control with r = 7/38, oscillations 
on the dip are significantly suppressed and the resolution 
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FIG. 2. (color online). Filter functions and coherence of the 
NV sensor for different controls, (a) Designed DD sequences, 
(b-e) Filter functions of DD controls with r = 1/3, 7/38, 5/18, 
and 3/10, respectively. The first dominant peaks are empha¬ 
sized in red. (f) Heights of the dominant peaks as functions 
of r. The dashed and solid curves correspond to the peaks 
at Lijtl{2n) = 15 and 5, respectively, (g-j) Coherence under 
DD controls with r = 1/3, r = 7/38, r = 5/18, and 3/10, 
respectively. The horizontal coordinate r donates t/n, where 
t is the total free evolution time. Note that the range of the 
horizontal axes in (h-j) is 1/3 of that in (g). The black circles 
represent experimental data. The red curves are calculations 
according to Eq. S16 multiplied a decay factor. 



FIG. 3. (color online). Calculations, (a) Coherence of the NV 
sensor with a nuclear spin nearby under DD controls with 
r = 1/3 (upper) and 3/10 (lower). The decoherence is ig¬ 
nored in the dashed curves, while it is taken into account in 
the solid curves, (b) Signal magnitudes of three nuclear spins 
resolved in experiment [Figs. 2(i,j)] as functions of r. (c-e) 
Signal magnitude maps of nuclear spins under DD controls 
with r = 1/3, 5/18, and 3/10, respectively. The chroma rep¬ 
resent the signal magnitudes, and the three color bars are the 
same. The dashed contours from the outside in correspond 
to signal magnitudes of 0.02, 0.1, and 0.5. The grey areas 
have the signal magnitudes below 0.01. (f) Signal magnitude 
difference between the cases of r = 5/18 and 3/10 by sub¬ 
tracting the former [Fig. 3(d)) from the latter [Fig. 3(e)[. The 
three dots represent three nuclear spins resolved in experi¬ 
ment [Figs. 2(i,j)]. 


is thus improved [Fig. 2(h)]. This allows for extracting 
the coupling information of single nuclear spins straight¬ 
forwardly from each dips. 

To selectively amplify the signal in III, the designed 
DD control with r = 5/18 is applied. The dips in /’ 
also arise, and they originate from the same nuclear spin 
as that in I. Under the designed DD control with r = 
3/10, a new coherence dip in IV emerges [Fig. 2(j)]. This 
dip, not resolved under CPMG control, is due to two 
remoter nuclear spins [31]. To sum up, the signals of 
nuclear spins for r 1/3 [Figs. 2(h-j)] move forward to 
1/3 positions in the time domain relative to the case of 
r = 1/3 [Fig. 2(g)], while the coherence times are not 
shortened as much, especially for r = 5/18 and 3/10. 
This relaxes the decoherence restriction and makes more 
and remoter nuclear spins detectable. 

To show the advantage of well-designed DD controls 


more clearly, we take a nuclear spin with d = Inm and 
6 = 60° for instance, where d and 9 represent the NV-^^C 
distance and the inclination angle of the NV-^^C vector 
with respect to the magnetic field B [Fig. 1(a)]. The 
calculations show clearly that the nuclear spin can be re¬ 
solved by the designed DD control with r = 3/10, but 
cannot by the original CPMG control (i.e., r = 1/3) un¬ 
der the same magnetic field [Fig. 3(a)]. In the following 
we elaborate the analysis by more calculations based on 
Eq. (1). The r dependence of the depths of coherence 
dips (regarded as signal magnitudes) for the three de¬ 
tected nuclear spins [Figs. 2(i,j)] is plotted in Fig. 3(b) 
[31]. It shows that optimal values of r vary among nuclear 
spins. Under the assumption that only dipole-dipole cou¬ 
plings are present between the sensor and nuclei, and that 
the dipole moment of the sensor is point-like, the signal 
magnitude of a nuclear spin in any spacial position can 
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be determined. In this way, signal magnitude maps are 
obtained. The maps plotted in Figs. 3(c-e) correspond 
to r = 1/3,5/18, and 3/10, respectively. The coordi¬ 
nates represent the longitudinal distance dy = d cos 9 and 
transverse distance d±_ = ds'md between the NV center 
and nuclei. The three maps indicate that designed DD 
controls with r = 5/18 and 3/10 surpass the undesigned 
one with r = 1/3 in detecting remoter nuclear spins, 
and this is consistent with the experimental results in 
Figs. 2(g,i,j). The discrimination between Figs. 3(d) and 
Figs. 3(e) is shown in Fig. 3(f). The DD control with 
r = 3/10 is superior in the positive areas. This dif¬ 
ference map implies that remoter nuclear spins can be 
observed with r = 3/10 rather than r = 5/18, and is 
consistent with the contrast between Figs. 2(i) and 2(j). 
It is evident that nuclear spins in different regions are 
selectively resolved through the adjusting process. The 
distances between the relevant nuclear spins are of the 
magnitude of several angstroms as schematically illus¬ 
trated in Fig. 1(a) [31], so the field of detection can be 
tuned on sub-nanoscale. This behavior may enable to¬ 
mography of the nuclear spin environment. 

Furthermore, this method has variations. We can 
group every five pulses together as a unit and modify 
the sequences in a similar fashion [Fig. 4(a)]. The filter 
function [Fig. 4(b)] is 




2 5ujt 5pujt 

Q —-cos- 

4t7, n 


5aujt \ ^ 


where n is divisible by 10 and 0 < p < q < 0.5. The 
dominant peaks of this function center at u}t/(2'K) = 
(2fc - l)n/10 (A: = 1,2, 3, • • •) for p ^ 1/5 or q ^ 2/5. 
The height of the fc-th dominant peak is 



{2.5-5q)T {2.5-5q)T 



FIG. 4. (color online). Designed DD control with five pulses 
per unit, (a) Pulse sequence, (b) Filter functions of DD 
controls with p = 3/20 and q = 17/40. The first dominant 
peak is emphasized in red. (c) Coherence under DD controls 
with p = 3/20 and q = 17/40. The horizontal coordinate r 
donates t/n, where t is the total free evolution time. The 
black circles represent experimental data. The red curves are 
calculations according to Eq. S16 multiplied a decay factor. 


to be resolved. This provides a route towards tomog¬ 
raphy of nuclear spin environment, and is inspiring in 
nanoscale magnetometry. It also extends the capabil¬ 
ity of using nuclear spins as quantum registers in quan¬ 
tum information processing. Additionally, the princi¬ 
ple of DD design may be applicable in other systems, 
such as trapped ion and superconducting qubit, for high- 
resolution noise spectroscopy. 


htiP, q) = ~ ^ [(2^ — l)7rp] + 2 cos [{2k — l)7rg]}^. 

( 5 ) 

The signals of nuclear spins for r 1/3 [Figs. 2(h-j)] 
move forward to 1/3 positions in the time domain rela¬ 
tive to the case of r = 1/3 [Fig. 2(g)] The signals move 
forward to 1/5 positions of that in Fig. 2(g) in the time 
domain. Using well-designed five-piece DD controls, the 
unwanted signal in P of Fig. 2(j) can be eliminated, and 
the signals in III and IV stand out [Fig. 4(c)]. 
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IV. CONCLUSION 

In conclusion, we report the first experimental demon¬ 
stration of DD design for identifying weakly coupled nu¬ 
clear spins in a solid system. We show that properly de¬ 
signed DD controls outperform the conventionally used 
CPMG controls in resolving weakly coupled individual 
nuclear spins, resulting in more and remoter nuclear spins 
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